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Abstract 

We present a study of s-channel single top quark production at the upgraded Tevatron pp collider, 
including the next-to-leading order (NLO) QCD corrections to the production and decay of the top 
quark. The "modified" narrow width approximation was adopted to preserve the spin of the top 
quark in its production and decay. We discuss the effect of the different 0{as) contributions on the 
inclusive cross section as well as various kinematical distributions after imposing the relevant cuts 
to select s-channel single top signal events. In particular the 0{as) decay contribution, while small 
in size, has a significant impact on several distributions. With the help of the best-jet algorithm 
to reconstruct the top quark we demonstrate that it is possible to study kinematical and spin 
correlations in s-channel single top events. We furthermore compare top quark spin measurements 
in two different basis and show how NLO corrections have to be taken into consideration in searches 
for the Higgs boson through W^H associated production at the Tevatron. 
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I. INTRODUCTION 



At the Fermilab Tevatron, the dominant single top quark production mechanisms are 



nant smgie to 

mm 



and the t-channel process bq — » tq ' 



mel process oq — » tg 



the s-channel process qq' W* — > tb 

(including bq' tq, also referred to as the W^-gluon fusion process) 
It has been shown that it is possible to disentangle the single top signals from various 
backgrounds^^. While top quark pair production is primarily a QCD process, single top 
quark production occurs dominantly in electroweak processes. It thus provides a probe of 
electroweak properties of the top quark. For example, one can study the Cabibbo-Kobayashi- 



Maskawa (CKM) matrix element 

(SM) Q,HQ,0,Q,Q,Hl2ii 



Ah, new W-t-b couplings beyond the Standard Model 
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25|, and CP violation 



H 22, Q 



in single 



top quark production and decay. The angular correlations among the deca^ 
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products of 
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polarized top quarks also provide a useful handle on these couplings j2 

0,0]. The t-channel production process can also provide us with information about the 
b quark distribution function inside the proton. The s-channel process is also a significant 



background to f 
decay H ^ bb 



igg s searches at the Tevatron in the production process qq' — ^ WH with 
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and other new physics search 



The expected single toT 

^ tw- 
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production rate, associated with a py-boson, via the process bg 
is too small to be observed at the Tevatron and therefore not considered in this paper. 

Searches for single top production were performed by both D0 Q and CDF in Run I. 
At the 95% confidence level, the D0 limit on the s-channel production is 17 picobarn (pb) 
and the CDF limit is 18 pb. At the same confidence level, the D0 limit on the t-channel 
production cross section is 22 pb and the CDF limit is 13 pb. The s-channel and t-channel 
single top quark processes can be probed separately at the Tevatron by taking advantage 
of 6-tagging using displaced vertices and differences in the kinematical distributions of the 
6-tagged jet. Usually, only one 6-tagged jet can be expected in the t-channel case while two 
6-tagged jets can be expected in the s-channel case. This is because the b quark produced 
with the top quark tends to be collinear with the initial gluon in the t-channel, therefore 
having a transverse momentum [px) that is too low to be observed. We shall focus on 
the s-channel single top quark process in this paper, the t-channel results will be discussed 
elsewhere 

The extraction of a signal is more challenging for single top production than for top pair 
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production since there are fewer objects in the final state and the overall event properties are 
less distinct from the large ly+jets background. Therefore, an accurate calculation including 
higher order QCD corrections is needed. The next-to-leading order (NLO)^(ct^) corrections 
to single top quark production have already been carried out in Refs. 
studies of the 0{as) corrections to the top quark decay have been done 
effect of the decay of the top quark and the top quark width have not been included in 
the previous studies of kinematics which focused on the differential cross section in on-shell 
production 



31 |lO| . S imilarly, 
Jl. However, the 



50l Isij ]. Since top quark production and decay do not occur in isolation from 
each other, a fully differential calculation should include both kinds of corrections together 

n 

and maintain all spin and angular correlations. In Ref. a formalism is presented which 
calculates the complete next-to-leading order QCD corrections to production and decay of 
single top quarks. In that calculation, the "modified" narrow width approximation is used 
to separate the production of the top quarks from their decay. Top quark width effects are 
also included. In order to study the spin properties of the top quark, the helicity amplitude 
method is used in the matrix element calculation, thus keeping full spin information for all 
final state objects. Ref. [3| adopts the phase space slicing method (with one cutoff scale) 
to calculate the NLO differential cross section. A preliminary phenomenology study using 
this method was presented in Ref. js^. A similar study using the subtraction method has 
also been published recently 5^ . 

In this study, we construct a next-to-leading order Monte Carlo calculation which treats 
0{as) corrections in a consistent way to both production and decay of the top quarks at 
the upgraded Fermilab Tevatron.* We examine various kinematical event properties in s- 
channel single top quark production. Our approach will allow experimentalists to compare 
their results with the theoretical predictions directly. A detailed phenomenological study, 
including background processes will be discussed elsewhere The paper is organized as 
follows: In Sec. EJ we present the inclusive cross section of the s-channel single top quark 
processes, including an examination of the theoretical cutoff dependence. The acceptance 
for s-channel single top events under various scenarios of kinematical cuts is presented in 
Sec. mil We study final state kinematical distributions in Sec. IIVI In Sec. |3 we consider 
the s-channel single top quark processes as a background to SM Higgs boson searches in the 



A similar study of top quark pair production and decay in e+e annihilation was given in Ref. [ssf. 



channel qq' — WH with decay H bb. Finally, we give our conclusions in Sec. IVII 

II. CROSS SECTION (INCLUSIVE RATE) 

In order to be able to compare theoretical predictions with experimental results it is 
important to not only determine the total production rate of single top events but also 
the rate of events passing kinematical cuts due to detector acceptances or the need for 
background suppression. Furthermore, reconstructing the top quark is critical to many of 
the physics goals of the Tevatron and the LHC. At the Tevatron Run II, experiments may 
hope for an accuracy of 2 GeV [3| in the top mass measurement. But the ability to achieve 
this accuracy depends on how well systematic effects - especially those associated with 
gluon radiation - are understood and controlled. It is therefore crucial to properly simulate 
final state particle distributions in single top events. In this section, we discuss inclusive 
production rates while the next section focuses on the acceptance of s-channel single top 
events under various kinematical cuts. The following section is devoted to a discussion of 
distributions of final state particles with an emphasis on the effects of gluon radiation. 

We present numerical results for the production of single top events considering the 
leptonic decay of the H^-boson from the top quark decay at the upgraded Tevatron (a 
1.96 TeV pp collider). Unless otherwise specified, we use the NLO parton distribution 
function set CTEQ6M j^, defined in the MS scheme, and the NLO (2-loop) running 
coupling as with A^ provided by the PDFs. For the CTEQ6M PDFs, A^ = 0.326 GeV 
for four active quark flavors. The values of the relevant electrowe ak p arameters are: a = 



1/137.0359895, = 1.16637 x 10~^GeV~'^, rrit = 178 GeV |58|, |59|, mw = 80.33 GeV, 
niz = 91.1867 GeV, and sin^ 9w = 0.231. Thus, the square of the weak gauge coupling is 
= A\/2m^G ^. Here, we focus our attention on the positively charged electron lepton 
(i.e., positron) only, though our analysis procedure also applies to the n lepton. Including 
the 0{as) corrections to qq\ the decay branching ratio of the W boson into leptons is 



Br(W — s> l^u) = 0.108 j60|. Unless otherwise specified, we will choose the top quark mass 
to be 178 GeV and the renormalization scale (/ir) as well as the factorization scale (fip) to 
be equal to the top quark mass. The top quark mass and scale dependences of single top 
events are investigated in the second and third part of this section. 



A. Theoretical Cutoff" Dependence 



The total cross section up to NLO can be obtained by adding up the contributions from 
the Born level contribution and the virtual correction and the real correction with either soft 
or hard gluons. We adopt the phase space slicing method with one cutoff scale (Smin) and 



the universal crossing functions 
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63[ to regularize both soft and coUinear singularities. 



The cutoff Smin serves to separate the real correction phase space into two regions: (1) the 
resolved region in which the amplitude has no divergences and can be integrated numerically 
by Monte Carlo (MC) methods and (2) the unresolved region in which the amplitude contains 
all the soft and collinear divergences and can be integrated out analytically. The resulting 
divergences will be canceled by the virtual corrections or absorbed into the parton structure 
functions in the factorization procedure. Since this cutoff is introduced in the calculation 
only for this technical reason and is unrelated to any physical quantity, the result must not 
depend on it. In other words, the sum of all contributions, virtual, resolved, and unresolved 
corrections must be independent of Smin- This is the case as long as Smin is small enough so 
that the soft and collinear approximations are valid. However, the numerical cancellation 
in the MC integration becomes unstable if Smm is too small. Furthermore, the jet-finding 
algorithm and other infrared-safe experimental observables also should be defined in a way 
that is consistent with the choice of Smin- In practice, one wants to choose the largest Smin 
possible within these constraints in order to minimize the running time of the MC integration 
program. For our study, we found a value of Smin = 5 GeV^ to be ap prop riate. A detailed 
discussion of the phase space slicing method can be found in the Ref. (52| |. 

That the total rate is indeed insensitive to the value of Smm for a large range is illustrated 
in Fig. which shows the sum of the virtual and unresolved real corrections {s + v) as well 
as the resolved contribution (real) to the s-channel single top quark process as a function of 
Smin- Although the contributions from the individual pieces vary, their sum (total) remains 
essentially constant for a large range of Smin- With the choices of yU/j, = yU^? = at the 
Tevatron, we obtain an inclusive cross section of the s-channel single top [t only) processes 
(with IV-boson decay branching ratio) as 47.9 fb, which agrees with Ref. j54| with the 
understanding that different electroweak parameters are used. 
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Figure 1: The theoretical cutoff Smin dependence of the inclusive s-channel single top quark cross 
section at the Tevatron with fin = ^ip = for = 178 GeV. The decay branching ratio 
t bW^{-^ e^u) has been included. 

B. Top Quark Mass and Renormalization/Factorization Scale Dependence 

To test the standard model and measure the CKM matrix element Vtb, one needs an 
accurate prediction of the single top quark production and decay to reduce the theoreti- 
cal uncertainty. Examining the top quark mass dependence can provide information about 
how accurately the top quark mass must be measured in tt events in order to reduce the 
uncertainty on the extraction of the CKM element Vtb from the measured single top cross 
section. Besides the top quark mass, the choices of renormalization and factorization scales 
also contribute to the uncertainty of the theoretical prediction. The renormalization scale 

is introduced when redefining the bare parameters in terms of the renormalized parame- 
ters, while the factorization scale fip is introduced when absorbing the collinear divergence 
into the parton distribution functions. Therefore, both /i/j and fip are unphysical and the 
final predictions should not depend on them. However, since we work at a fixed order in 
perturbation theory, we indeed see dependences of the predicted cross section on fiR and 
fip- The change due to varying the scale is formally of higher order. Since the single top 
rate is small at the Tevatron, it is very important to reduce the scale uncertainty in order 
to compare the theory prediction with experimental data. Here, we examine the top quark 
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Figure 2: Top quark mass as well as renormalization and factorization scale dependence of the 
inclusive s-channel single top quark cross section at Tevatron. The decay branching ratio of t ^ 
bW^{-^ e'^v) has been included. 

mass and the scale dependence of s-channel single top quark events. 

As shown in Fig. |2l the cross section changes by about ±10% when the top quark mass 
mt is varied by its current uncertainty of about =f5 GeV. Measuring the top quark mass to 
an uncertainty of 1 — 2 GeV will thus reduce the theoretical uncertainty on the single top 
cross section. The reduction in uncertainty of the s-channel single top production rate will 
improve the measurement of the CKM matrix element Vtb- 

To examine the scale dependence of the s-channel single top quark production rate, we 
show in Fig.Elthe results of two typical scales; one is the top quark mass (/ii? = = rrit), 
shown as the solid-line, and another is the total invariant mass of the event {fip = = 
V^), shown as the dashed-line. For the decay of top quark, we take fij^"^ = nit which 
gives similar results as the choice of ^'^'^ = My/- The band constrained by these two 
scales represents a range of uncertainty due to the NLO predictions. The usual practice 
for estimating the yet-to-be calculated higher order QCD correction to a perturbative cross 
section is to vary around the typical scale by a factor of 2, though in principle the "best" 
scale to be used for estimation cannot be determined without completing the higher order 
calculation. In Fig. IHlwe show the total cross section of the s-channel single top production 
for a range of the scale /i^ which, for simplicity, is set to be equal to ^r. We have multiplied 
a constant factor (the "K-factor") of 1.52 to the LO cross sections (shown as the dashed-line) 
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Figure 3: Inclusive s-channel single top quark cross section at Tevatron for mj = 178 GeV, versus 
the ratio of the factorization scale fip to its typical value /xq, where /io = mt and /zq = 
respectively. The decay branching ratio of t — > bW~^{^ e'^v) has been included. 

in order to compare to the NLO ones (shown as the solid-hne). It is clear that the NLO 
calculation reduces the scale dependence. For example, around = mt, the LO rate varies 
by +8.9% to —5.7% when scale is changed from fipf^ to 2fiF, while the NLO rate varies by 
+7.6% to —4.7%. Similar results also hold for varying the scale around fip = V^- 

C. Inclusive Cross Section 

Due to the electroweak nature of single top quark production, the higher-order QCD 
corrections can be divided into three separate gauge invariant sets: the corrections from 
the initial state of top quark production (INIT), the corrections from the final state of 
top quark and b quark (FINAL), and the corrections from the top quark decay (SPEC). 
The explicit diagrams and definitions for these three sets can be found in Ref. j52|. The 
inclusive cross section as well as the individual 0{as) contributions are listed in Table IB 
The effect of the finite widths of the top quark and l^-boson has been included. We use 
a "modified" narrow width approximation (MNWA) in the calculation instead of the usual 
narrow width approximation (NWA). In the usual NWA, the effect of the Breit-Wigner 
resonance propagator reduces to a delta-function in the limit of vanishing top quark width. 
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I.e. 




(p2 — m^)^ + m^Fj 



1 



6 {p^ - ml) . 



(1) 



Therefore, the invariant mass of the top quark decay particles will be exactly equal to 
rrit (a fixed value). To model the reconstructed top quark invariant mass from its decay 
particles with a Breit-Wigner resonance shape to reflect the non-vanishing decay width of 
the top quark (for being an unstable resonance), we introduced the MNWA method in our 
numerical calculation in which we generate a Breit-Wigner distribution for the top quark 
invariant mass in the phase space generator and then calculate the squared matrix element 
with rrit being the invariant mass generated by the phase space generator on the event-by- 
event basis. We find that the total event rate and the distributions of various kinematics 
variables (except the distribution of the reconstructed top quark invariant mass) calculated 
using the "modified NWA" method agree well with that using the NWA method. In the 
NWA method, the reconstructed top quark invariant mass distribution is a delta-function, 
while in the "modified NWA" method, it is almost a Breit-Wigner distribution. The reason 
that the "modified NWA" method does not generate a perfect Breit-Wigner shape in the 
distribution of the top quark invariant mass is because the initial state parton luminosities 
(predominantly due to valence quarks) for the s-channel single-top process drop rapidly 
at the relevant Bjorken-x range, where (x) ^ ^ ~ 0.1. As shown in Table HI the total 
0{as) QCD corrections increase the Born level cross section of the single top quark event 
by 54%. The INIT correction dominates over the FINAL and SDEC corrections, due to 
the enhancement from coUinear physics. The contribution from the SDEC corrections is 
furthermore suppressed by almost a factor of 20 compared to the INIT corrections. This 
can be understood from Eqs. (146)-(149) in Ref. js^ and is explained below. 

When calculating the NLO QCD corrections to the top quark decay process, we must use 
the NLO top quark decay width in order to restore the correct top quark decay branching 
ratio. With the narrow width approximation, the NLO differential cross section of the top 
quark decay process is given by 





TT 



(2) 



X 



where 




Here r° denotes the Born level top quark decay width, and T] denotes the 0{as) correction to 
the top quark decay width. In our calculation, = 1.695 GeV and rf^*-* = 1.558 GeV 
The 0{as) SDEC correction includes only the 0{as) QCD correction in Eq. (0), 

d^ll^l") = dal-'^-'- X ^ X dTt-'. (3) 



We obtain 



t 



IT TT TT TX f _ T\ 



0{al). (4) 



Expanding Eq. ^ up to 0(as), we get 

X ^ X ^ X c/r;^^^"^ + daP^o°^"^*'°" X ^ X drf + 0(a2). (5) 

The first term in Eq. belongs to the Born level contribution, and both the second and 
third terms belong to the 0(0:5) decay contribution. The net contribution of the second 
and third terms should be exactly zero in order to respect the top quark decay branching 
ratio (to be one), after we integrate out the phase space of the top quark decay products. 
As shown in Table HI there is a small remaining SDEC contribution, which is due to higher 
order contributions: 

• The 0{a1) term, which is dropped in the expansion in Eq. (0), will give about a 1% 
contribution. 

• A mismatch exists between J dV'^'^^^ and F^, since we use a Breit-Wigner distribution 
of the top quark mass modulated by its decay width in the MNWA method rather 
than a fixed mass (which indicates an on-shell top quark). 

This SDEC contribution will become sizable once we use a jet finding algorithm to define 
infrared-safe observables and impose the kinematical cuts. This is especially the case for 
three-jet events which can only come from the real emission corrections. However, the decay 
contribution cannot be very large because the large top quark mass has already set the scale 
for the decay of the top quark which will be shown later. 
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uross bection (ibj 


Fraction of NLO (%) 


Born Level 


31.2 


65.0 


INIT 


10.7 


22.3 


FINAL 


5.5 


11.5 


SDEC 


0.57 


1.19 


0{as) sum 


16.8 


35.0 


NLO 


47.9 


100 



Table I: Inclusive single top production cross section for different subprocesses, including the top 
quark decay branching ratio t bW^{-^ e'^v). 

III. SINGLE TOP ACCEPTANCE STUDIES 



In this section we explore the final state objects of s-channel single top events. Although 
the ly-boson from the top quark decay can decay in both leptonic and hadronic modes, we 
only focus on the leptonic decay mode in this study because the all-jet production mode of 
single top events is difficult to observe experimentally. Therefore, the signature of s-channel 
single top events that is accessible experimentally consists of one charged lepton, missing 
transverse energy, and two or three jets. Since we are studying the effect of NLO QCD 
radiative corrections on the production rate and the kinematical distributions of single top 
events at the parton level in this paper, we are not including any detector effects, such as jet 
energy resolution or b-tagging efficiency. Only the kinematical acceptance of the detector is 
considered. Since the Tevatron is a pp collider, and pp is a CP-even state, the production 
cross sections for hi[j) at the Tevatron are the same as those for th{j) (when ignoring the 
small CP violation effect induced by the CKM mixing matrix in quarks). For simplicity we 
therefore only show distributions of s-channel single top events in this work in which a t 
quark (not including t) decays into W^{-^ l^v) and h quark. 

To meaningfully discuss the effect of gluon radiation in single top events, we have to 
define a jet as an infrared-safe observable. The numerical stab ility of several jet algorithms 
commonly used in experiments has been investigated in Ref. In this study, we adopt 
the cone-jet algorithm 6^ which is organized as follows: 

1. Build a sorted list of all clusters (in this case the final state h and h quarks plus a light 
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quark or gluon from the 0{as) correction), decreasing in transverse energy Enp. 

2. Select the highest cluster from the cluster list and draw a cone of radius R around 
its axis. Call this a jet and calculate the transverse jet energy and a new jet axis by 
summing over the 4-momenta of all clusters inside the cone. 

3. Remove all clusters in the cone from the cluster list and move the newly constructed 
jet to the jet list. 

4. Apply the appropriate minimum transverse energy and rapidity cuts to the entries in 
the jet list to create the final state hst of jets. 

More specifically, we adopt the i?-scheme cone-jet (4-momenta of particles in one cone are 
simply added to form a jet) with radius R = y^Arp~+~A0^ to define 6, b and possibly extra 
g or q (or q) jets, where Arj and A(f) are the separation of particles in the pseudo-rapidity rj 
and the azimuthal angle 0. For reference, we shall consider both R = 0.5 or R = 1.0. The 
same -R-separation will also be applied to the separation between the lepton and each jet. 
The kinematical cuts imposed on the final state objects are: 

> 15GeV , \r]i\ < t]™'^^ 

> 15 GeV , 

> 15 GeV , \r]j\ < r/f'^^ 

ARij > Rcut , '^Rjj > Rcut (6) 

where the jet cuts are applied to both the b- and fe-jet as well as any gluon or light anti-quark 
jet in the final state. Two lepton pseudo-rapidity cuts are considered here: a loose version 
with rj^"-^ = 2.5; and a tight version with 7^™'^^ = 1.0. Similarly, loose and tight cuts are 
also considered for the jet pseudo-rapidity, rjj^'^^ = 3.0 and rij^"-^ = 2.0, respectively. The 
minimum transverse energy cuts on the lepton as well as the jets is 15 GeV. Each event is 
furthermore required to have at least one lepton and two jets passing all selection criteria. 
The cut on the separation in R between the lepton and the jets as well as between different 
jets is given by Rcut which is is chosen to be 0.5 (or 1.0). 

In Table ITTl we show the single top production cross sections in femtobarns (fb) as well as 
acceptances for different subprocesses for several sets of cuts. We apply the Et cuts listed 
in Eq. © and study three separate sets of values: 
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Cross section (fb) 


Acceptance (%) 






LO 


NLO 


INIT 


FINAL 


SDEC 


LO 


NLO 


(a) 


tb + tbj 


22.7 


32.3 


6.6 


3.0 


0.16 


72 


67.4 




t6 + tbj, Etj > 30 GeV 


17.8 


24.8 


5.5 


1.7 


-0.16 


57 


52 




tbj 




10.6 


6.2 


2.6 


1.8 




22.1 




tbj, Exj > 30 GeV 




3.5 


2.2 


0.92 


0.33 




7.3 


(b) 


tb + t6j 


19.0 


21.7 


2.0 


1.4 


-0.54 


60.9 


45.3 




tb + tbj, Etj > 30 GeV 


14.8 


16.7 


1.8 


0.70 


-0.48 


47.4 


34.9 




tbj 




5.6 


3.9 


1.1 


0.56 




11.7 




tbj, Etj > 30 GeV 




1.9 


1.4 


0.40 


0.08 




4.0 


(c) 


tb + tbj 


14.7 


21.4 


4.4 


2.2 


0.13 


47 


44.7 




tb + tbj, Etj > 30 GeV 


12.2 


16.9 


3.6 


1.3 


-0.17 


39 


35.3 




tbj 




6.4 


3.6 


1.7 


1.1 




13.4 




tbj, Erj > 30 GeV 




2.3 


1.5 


0.65 


0.23 




4.8 



Table II: The s-channel single top production cross section (in fb) and acceptance at the Tevatron for 
different subprocesses under various scenario: (a) is for the loose cuts with rj'^"-^ = 2.5, r/™"^ = 3.0, 
and Rcut = 0.5, (b) is also for the loose cuts, but with a larger jet clustering cone size (and jet- 
lepton separation cut) of Rcut = 1-0, and (c) is for the tight cuts with r/™"^ = 1.0, ■q'j^"'^ = 2.0, and 
Rcut = 0.5. The first five columns show the cross section at Born level and NLO as well as the 
different 0{as) contributions. The last two columns show the acceptance at Born level and NLO. 
The decay branching ratio t bW~^{-^ e^v) is included. 

(a) loose cuts with small Rcut- '<lT°'^ = 2.5, rfj^"-^ = 3.0, and Rent = 0.5, 

(b) loose cuts with large Rcut'- vr"'^ = 2.5, rjj'"-^ = 3.0, and Rcut = 1-0, 

(c) tight cuts with small Rcut- vr""" = l-O. VT""" = 2-0, and Rent = 0.5. 

A larger value for Rcut reduces the acceptance significantly mainly because more events fail 
the lepton-jet separation cut. While this is only a 16% reduction at Born level, the difference 
grows to 33% at NLO. Hence, a smaller cone size is preferred to keep the acceptance at a 
high level. For events with at least three jets, imposing a harder cut on the transverse 
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momentum of the third jet {Etj > 30 GeV) only decreases the contribution from INIT and 
FINAL corrections by a factor of 3. By contrast, this cut reduced the contribution from the 
SDEC correction by a factor of more than five because the top quark mass sets the scale for 
the top quark decay contribution rather than the invariant mass of system. A tighter cut 
on jet Et also reduces the relative SDEC contribution in the 3-jet bin. The Et spectrum of 
the gluon from the SDEC contribution is softer because the available phase space is limited 
by the top quark mass, which can also be seen below in Fig. E| 

We note that the acceptances at Born level and NLO for inclusive 2-jet events differ 
by 7% for the loose set of cuts (case (a) in Table |n} and by 5% for the tight set of cuts 
(case (c) in Table ITT|l . In both cases, the acceptance is reduced at NLO compared to Born 
level. This reduction grows to over 30% when a larger value for Rcut of 1.0 is used. Since 
the acceptances are quite different between Born level and NLO, we cannot use a constant 
K-factor to scale the Born level distributions to NLO. The kinematical differences between 
Born level and NLO that result in different acceptances will need to be taken into account 
properly in the measurement of Vtb- 

Fig. m shows how the observed cross section changes as a function of the jet Et cut when 
applying the loose set of cuts, including a requirement of at least 2 jets. The cross section 
(and thus the acceptance) doesn't change very much until the jet E^ threshold reaches about 
25 GeV because the b- and 6-jets typically have high Et- The figure also shows that a jet 
pseudo-rapidity cut of \ri\ < 2 or |?7| < 3 does not impact the overall acceptance; the cross 
section only decreases significantly when restricting jets to the very central region (|?7| < 1). 
This implies not only that t] cuts on jets should be large when separating single top events 
from the large backgrounds, but also that 6-tagging needs to be efficient over a large rj range 
(at least \ri\ < 2) to be able to detect s-channel single top signal events. 

The dependence of the fraction of 2-jet events and 3-jet events on the jet Et cut is also 
shown in Fig. EJ At Born level, there are only 2-jet events, whereas 0{as) corrections can 
produce an additional soft jet. The fraction of events with these additional jets is low only 
for very high jet Et thresholds. For typical jet Et thresholds considered by experiments of 
15 GeV to 25 GeV these jets add a significant contribution. As expected, the effect is not 
quite as large when only jets within a very small rj range are considered because the extra 
jet typically has higher r]. In order to study 0{as) effects it is thus important to set the jet 
T] cut as high as possible and the jet Et cut as low as possible. 
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Figure 4: Cross section and fraction of 3-jet events at NLO for varying jet Et cuts, requiring only 
that rijets ^ 2, and not making any cuts on the electron or neutrino. The left-hand figure shows 
the total cross section for events with 2 or 3 jets as a function of the jet Et cut for three different 
jet pseudo-rapidity cuts. The right-hand figure shows the fraction of 3-jet events as a function of 
the jet Et for different jet pseudo-rapidity cuts. The lowest threshold considered is 15 GeV. 

As mentioned before, the event rate for single top events is small and it is important 
for experiments to maximize their acceptance. We will thus use the loose set of cut values 
for the following discussion: 7^™'^^ = 2.5, rjj^""^ = 3.0, and Rcut = 0.5, EJPj^ = 15 GeV, cf. 
Eq. ©. 

IV. SINGLE TOP EVENT DISTRIBUTIONS 

In this section, we study the kinematical properties of single top quark events. For the 
s-channel process at the Born level, two 6-jets (the b from the top quark decay and the b 
produced with the top) appear in the final state and cannot be distinguished from each other 
experimentally. A prescription is needed to identify which of the two jets corresponds to the 
b quark from the top decay. At NLO, an addition jet is radiated, which further complicates 
the reconstruction of the top quark. This is because the additional jet can come from either 
the production or the decay of the top quark. Production-stage emission occurs before the 
top quark goes on shell and decay-stage emission occurs only after the top quark goes on 
shell. In production emission events, the W boson and b quark momenta will combine to 
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Figure 5: The transverse momentum pT (left) and pseudo-rapidity r] (middle) distributions of the 
electron and the missing transverse energy (right) after selection cuts, comparing Born level to 
O(q^) and full NLO. 

give the top quark momentum, while in decay emission events the gluon momentum must 
also be include to reconstruct the top quark momentum. To find the best prescription of 
classifying the gluon jet correctly, we first examine various kinematical distributions of the 
final state particles. We then investigate two top quark reconstruction prescriptions: the 
leading jet algorithm and the best-jet algorithm. Having chosen a prescription, the effects 
of NLO corrections on distributions concerning the reconstructed top quark are studied, in 
particular spin correlations between the final state particles. Finally, we explore the impact 
of the radiated jet in exclusive 3-jet events. 



A. Final State Object Distributions 

In this section we examine various kinematical distributions of final state objects after 
event reconstruction and after applying the loose set of cuts, cf. Table ID^a) and Eq. 0. 
We study inclusive 2-jet events in this section because they give more reliable infrared-safe 
predictions. 

Fig. El shows the transverse momentum of the electron and the missing transverse mo- 
mentum As expected because they are leptons and not quarks, the change in shape 
when going from Born level to 0{as) is not very large for the electron or the neutrino where 
the 0{as) corrections have the same shape and peak position as the Born level contribution. 
The pseudo-rapidity distribution of the electron is also given in Fig. This distribution 





Jet, pj [GeV] Jetg [GeV] 

(c) (d) 

Figure 6: Transverse momentum px of the b jet (left) and the b jet (right) after selection cuts, 
comparing Born level to 0{as) corrections. The top row shows the NLO, Born, and sum of 0{as) 
contributions, the bottom row shows the individual 0{as) contributions. 

widens at 0{as) and shifts slightly to more positive values of 77, but again the effect is not 
very large. We note that the peak position of the distribution is at a higher value than 
that of the electron pj- because the neutrino from the VT-boson decay moves preferentially 
along the direction of the top quark. This is due to the left-handed nature of the charged 
current interaction and can easily be seen when examining the spin correlation between the 
charged lepton and the top quark in the top quark rest frame. We will comment more on 
this subject in Sec. IIVC'I 

Compared to the lepton and neutrino, the effect of the 0{as) corrections on the recon- 
structed 6-jet and 6-jet is more pronounced. Fig. IHl shows a comparison of the distribu- 
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tions of these two jets between the Born level and the 0{as) contributions. Though it is not 
possible to distinguish between the b- and 6-jet experimentally, it is nevertheless instructive 
to consider their Et distributions individually. At the Born level, the transverse momentum 
of the 6-jet peaks higher and drops off faster than that of the 6-jet, resulting in similar mean 
Et for the two (68 GeV for the 6-jet and 71 GeV for the fo-jet). Since the 6-jet is produced 
in association with the heavy top quark, it has a long tail into the high px region to balance 
the top quark. The px distribution of the 6-jet (from top quark decay) on the other hand is 
predominantly controlled by the top quark mass and therefore peaks at ~ mt/S. The NLO 
QCD corrections broaden the LO transverse momentum distributions and shift the peak 
positions to lower values. The result for the mean of the distribution depends on the cuts 
that are applied because the different 0{as) corrections have different effects; in the case 
of the loose cuts both distributions have a mean of 66 GeV at 0{as). The INIT correction 
shifts the mean of the 6-jet pt distribution up while the FINAL and SDEC corrections tend 
to shift it down. The 6-jet p^ distribution receives a large contribution from the FINAL 
correction and a negligible contribution from the SDEC correction, as expected. When a 
gluon is radiated from {t, h) quark line, it prefers to move along the 6-jet direction due the 
coUinear enhancement and therefore shifts the 6-jet px distribution to the small pt region, 
as shown in Fig. IHfd). For the same reason, the SDEC correction to the 6-jet pt distribution 
also peaks around 60 GeV (the Born level peak position), but it provides a negative contri- 
bution, cf. Fig. inH^c). Furthermore, because the FINAL correction is larger in size than the 
SDEC correction, the resulting shift to the low pt region is stronger for the 6-jet than the 
6-jet. 

This shift in the rj distribution of the 6-jet, when comparing the Born level to 0(as), is 
also evident in Fig. [71 where in particular the FINAL correction shifts the distribution to 
more central pseudo-rapidities. The top quark is so heavy that it is mostly produced in 
the central rapidity region and thus its partner 6-jet also peaks around a pseudo-rapidity 
of zero. We note that the pseudo-rapidity distribution of the fc-jet is slightly asymmetric at 
the Born level, denoted by the red dotted line in Fig. [71 favoring negative pseudo-rapidity. 
This is similar to the lepton pseudo-rapidity asymmetry observed in the pp — > — > l^u 
events due to the ratio of down-quark and up-quark parton distributions inside the proton 
and anti-proton. The 6-jet coming from the top quark decay follows the top quark moving 
direction and therefore favors a slightly positive pseudo-rapidity, denoted by the solid line 
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Figure 7: Pseudo-rapidity rj of the 6-jet (left) and the 6-jet (right) after selection cuts, comparing 
Born level to 0{as) corrections. 
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Figure 8: Total event transverse energy Ht after selection cuts, comparing Born level to 0{as) 
corrections. 

in the Fig. [7| The 0{as) corrections, in particular the FINAL correction, shift the &-jet 
to more central regions. The shape of the fe-jet pseudo-rapidity distribution remains almost 
unchanged compared to Born level because it comes from the top quark decay. Therefore we 
expect to see an asymmetry between r/^ and rji, shown below after full event reconstruction. 

The impact that the different 0{as) corrections have on the pt of the jets is also reflected 
in event-wide energy variables such as the total transverse energy (Ht) in the event, defined 
as 

if^=p^-f*-+^y + (7) 

jets 

The distribution of Ht for s-channel single top events is shown in Fig. |H1 Similar to the b- 
and 6-jet pj- distributions, the FINAL and SDEC contributions shift the Ht distributions 
down, while the INIT contribution shifts it up. This is expected because the INIT correction 
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contributes additional energy to the event in the form of a third jet. 



B. Event Reconstruction 

When selecting single top events we would like to take advantage not only of simple single- 
object kinematics but also of correlations between objects. For that we need to reconstruct 
the event completely, not just the final state jets but also intermediate particles, in particular 
the VT-boson and the top quark. 

The VT-boson can be reconstructed from the final state electron and the observed missing 
transverse energy ^t- The lack of information about the beam-direction component of the 
neutrino momentum (p^) that would prevent this reconstruction is overcome by requiring 
the invariant mass of the electron-neutrino system to be equal to the mass of the 14^-boson. 
This additional constraint yields two possible solutions for p^, and typically, both of them 
are physical solutions for a signal event. In our analysis, we follow the prescription given 
in Ref. jo^] to choose the solution which has the smaller This picks the correct in 
about 70% of the events passing the loose set of cuts. 

To reconstruct the top quark, the reconstructed H^-boson then needs to be combined with 
the 6-jet from the top quark decay. The challenge that has to be overcome in this case is 
to identify the correct jet as the 6-jet. Several different methods have been used in the past 
to select this jet. The simplest method is to choose the highest-E'T jet in the event (leading 
jet). Another possibility that can be used when 6-tagging is available experimentally is to 
choose the 6-tagged jet in the event. However, neither of these methods is very successful 
in s-channel single top production because both the h quark from the top decay and the h 
quark produced with the top quark have high Et and central pseudo-rapidities. While the 
Et distributions are different for h- and 6-jets, Fig. [Ul shows that they differ both in peak 
position and in shape. For jet Et values above 100 GeV, the leading jet is actually more 
likely to be from the h quark. The leading jet corresponds to the h from the top about 55% 
of the time in the sample after loose cuts. 

A more effective method to identify which jet corresponds to the h quark from the top 
decay is to make use of the known top quark mass that is measured in ti events. In this 
"best-jet" algorithm, the Wj combination that gives an invariant mass closest to the true 
top mass is chosen as the reconstructed top quark. The jet thus identified is called "best- 
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Figure 9: Transverse momentum distribution of the b jet from the top quark decay, for ah events 
that pass the loose set of cuts (sohd Une), only for those events in which the b jet is also the leading 
jet (dashed line), and only for events in which the b jet is the best-jet (dotted line). 




If an additional jet is produced in the decay 0{as) contribution, then the best-jet method 
will not be able to reconstruct the top quark correctly. We therefore need to extend the 
best-jet algorithm to also include 2-jet pairs when forming a top quark. In other words, 
candidates include not only Wj but also Wjj. The list of systems for which the invariant 
mass is evaluated and from which the best-jet is chosen thus consists of: {W,jeti), {W,jet2), 
(W,jeti,jet3), (W,jet2,jet-^).'' The jet (or 2-jet system) that is thus chosen will be called 
the "best-jet" . Once we have identified the fe-jet from the top quark decay, we can assign one 
of the other jets to the b quark produced with the top. We choose this "non-best- jet" also 
from the leading two jets and assign the label to whichever of the two jets was not identified 
as the "best-jet" . 

The best-jet algorithm identifies the 6-jet properly about 80% of the time. The effec- 
tiveness of this method is mainly limited by the efficiency of the VT-boson identification 
method; if the H^-boson is not reconstructed properly, then identifying the b quark from the 
top decay becomes a random pick. Fig. IHl shows the dependence of this efficiency on the 
transverse momentum of the 6-jet. As expected, the leading jet corresponds to the b quark 
from the top quark decay mostly when the transverse momentum of the 6-jet is very large. 
The best-jet algorithm in comparison shows high efficiency for all transverse momenta. 

^ The combination {W,jeti,jet2) is not considered because it is very unlikely for the gluon jet to be one of 
the leading two jets due to its generally low px and because the gluon jet will not be b-tagged. 
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Figure 10: Invariant mass of the reconstructed W and a jet object. This jet object is: the jet 
containing the b quark from the top decay, using parton information (sohd Hne); the best-jet 
identified using the method outhned in the main text (dashed hne); the leading jet (dotted hne). 

Fig. shows the invariant mass of the reconstructed top quark, comparing the different 
methods to identify the b quark from the top decay. For the 6-jet curve in Fig. parton- 
level information is used to identify which of the final state jets contains the b quark from 
the top decay. If this b quark produced a gluon in decay-stage emission, this resulting extra 
jet is also included when reconstructing the top quark. For the leading jet curve, only the 
highest Et jet in each event is used. 

The top quark width is larger than it would be at parton level even though no smearing 
was applied in this study. This is the result of using the reconstructed kinematics of the W- 
boson (in particular the neutrino z-momentum) rather than parton-level information. Since 
the same reconstructed W^-boson was used in each case, differences between the individual 
curves are solely due to which jet is chosen to reconstruct the top quark. Furthermore, 
because parton-level information is used for the 6-jet curve, it functions as a reference and 
upper limit for the other methods. 

As expected, using the leading jet gives a peak at the invariant mass of the top quark with 
a height of about half the fe-jet peak. The best-jet algorithm by contrast shows approximately 
the same width and height as the 6-jet curve. This is the result of two competing effects, 
to be discussed in order. On the one hand the best-jet algorithm identifies the correct jet 
only in some fraction of the events because the W^-boson itself is mis-reconstructed part of 
the time, thus reducing the height of the curve. On the other hand the algorithm chooses 
an invariant mass close to 178 GeV and thus tends to bring mis-reconstructed events closer 
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Figure 11: Transverse momentum and rapidity distributions of top quark at the parton level 
without any kinematical cuts (a, 6), and after reconstructing the top quark from the final state 
particles with kinematic cuts {c,d), comparing Born level to 0{as) corrections. 

to 178 GeV. The latter effect is slightly larger than the former, hence the best-jet curve is 
slightly higher than the 6-jet curve in Fig. ^1 

Fig. shows the transverse momentum and rapidity distributions of the top quark at 
the parton level without any kinematic cuts (a, b) and after reconstructing the top quark 
from the the final state particles {W boson and the best-jet) with kinematic cuts (c, d). Our 
parton level result, without imposing any kinematical cut, is in good agreement with that in 
Ref. After event reconstruction, the 0{as) QCD corrections change the parton level top 
quark transverse momentum distribution, especially in the low pt region. The corrections 
also shift the top quark rapidity distribution to the more central region. We note that using 
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Figure 12: Transverse momentum difference between the best-jet and the non-best-jet (left) and 
pseudo-rapidity difference between the two (right) after selection cuts, comparing Born level to 
0(qs) corrections. 

the best-jet algorithm to reconstruct top quark results in distributions that are very similar 
to those obtained using the true 6-jet from top quark decay. 

C. Kinematical and Spin Correlations 

Having identified the 6-jet from the top quark decay through the best-jet algorithm, we 
can now study correlations expected from event kinematics. In s-channel single top events, 
there is a strong correlation between the kinematics of the b and b quarks at the Born level. 
This correlation is modified by gluon radiation in the production and decay of the top quark. 
Fig. 1121 shows how 0{as) effects change the momentum difference and the pseudo-rapidity 
difference between the best-jet and the non-best-jet. 

The transverse momentum difference is not affected very much; while the 0{as) correc- 
tions affect both the 6-jet and the 6-jet and tend to make the distribution broader, this is 
a small change compared to the overall width of the distribution. The 0{as) corrections 
have a larger affect on the pseudo-rapidity difference, shifting the distribution to more cen- 
tral values and making it broader. Both effects are mainly due to the FINAL and SDEC 
corrections, which tend to weaken the correlation between the 6-jet and the 6-jet. 
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Figure 13: Illustration of the correlation between the charged lepton from the top decay and the 
top quark spin, in the top quark rest frame. 

1. Top Polarization 



In the SM, the top quark produced in single top events is highly polarized, and this 
polarization can in principle be measured. Two definitions for the polarization have been 
studied in the literature, differing by the reference frame used to define the polarization: 
One calculation uses the helicity basis, another the so-called "optimal" basis j68|, |69|. In the 
more common helicity basis the top quark spin is measured along the top quark direction of 
motion in the center of mass (cm.) frame which is chosen as the frame of the (reconstructed 
top quark, non-best-jet) system after event reconstruction. In the optimal basis (beamline 
basis) we can maximize the spin correlations by taking advantage of the fact that the top 
quark produced through the s-channel single top quark processes is almost 100% polarized 
along the direction of the d-tjpe quark, which comes predominantly from the anti-proton. 
In the discussion below, we will examine the polarization of single top quark events in both 
cases. The particular strengths and weaknesses of both methods were discussed in Ref. [2^ . 

It has been pointed out that among the decay products of top quark the charged lepton 
is maximally correlated with the top quark spin 0, 0|. We can thus obtain the most 
distinctive distribution by plotting the angle between the spin axis and the charged lepton 
in the top quark rest frame. Before examining this spin correlation effect for a particular 
process or in a particular basis, it is worthwhile to give a schematic picture. Due to helicity 
conservation and the left-handedness of charged current interactions, the VT-boson from the 
top quark decay can be either longitudinally or left-handed polarized. This is true for a 
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massless b quark and diagrammatically shown in Fig. ^| The figure shows the preferred 
moving direction of the lepton from a polarized ly-boson in the rest frame of a polarized top 
quark. When the is left-handed polarized, the fact that the b quark must be left-handed 
forces it to move along the direction of the top quark polarization, cf. the upper part of 
Fig. I13r a). The thus moves against this direction. When the decays, the charged 
anti- lepton (e"*") must be right-handed, hence it prefers to move against the direction, in 
the same direction as the top quark polarization. When the is longitudinally polarized, 
it prefers to move in the same direction as the top spin, cf. the lower part of Fig. IT^ a). Its 
decay products prefer to align along the polarization, and since the is boosted in 
the direction of the top quark polarization, the charged anti-lepton again prefers to move 
along the top quark spin axis. Since in both cases the charged anti-lepton moves against 
the top quark moving direction, the neutrino from the decay will be harder than the 
charged anti-lepton, as shown in Fig. El Similarly, the preferred moving direction of the 
charged anti-lepton from a right-handed polarized top quark is illustrated in Fig. IT^ b). For 
simplicity, we will continue to use the phrase "charged lepton" when referring to the charged 
anti-lepton from the W decay. 

In the helicity basis, the polarization of the top quark is examined as the angular dis- 
tribution (cos^/ie/) of the lepton in the cm. frame of the incoming partons relative to the 
moving direction of the top quark in the same frame. The angular correlation in this frame 
is given by 

cose,,i = ^;^, (8) 
\ptM\ 

where pt is the top quark three-momentum defined in the cm. frame of the two incoming 
partons, and is the charged lepton three-momentum defined in the rest frame of the top 
quark. For a left-handed top quark, the angular correlation of the lepton i'^ is given by 
(1 — cos 9hei)/ 2, and for a right-handed top quark, it is (1 + cos 9hei)/'^- Fig. IT^ shows that 
this linear relationship for cos 6hei is indeed a valid description for s-channel single top events 
at the parton level. The figure also shows that the top quark is not completely polarized 
in the helicity basis, and that this polarization is weakened further when including 0{as) 
corrections. 

The cm. frame after event reconstruction is chosen as the frame of the (reconstructed 
top quark, non-best-jet) system. After reconstruction the effect of the lepton-jet separation 
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Figure 14: Top quark polarization in the helicity basis using the full parton information (left) and 
after event reconstruction (right) with selection cuts, comparing Born level to 0{as) corrections. 

cut can be seen as a drop-off of the cos 6 hei distribution close to a value of —1. 

In the "optimal" basis, the relevant angular correlation for the s-channel process is cos 9opt, 
defined as 



(9) 



COS 9opt 

where is the anti-proton three-momentum in the top quark rest frame and is the lepton 
three-momentum in the top quark rest frame. In this analysis, we orient the coordinate 
system such that the protons travel in the positive z direction; the antiprotons travel in the 
negative z direction. For a top quark polarized along the anti-proton moving direction, the 
angular distribution of the lepton is (1 + cos 9opt)/2, while for a top quark polarized along 
the proton moving direction it is (1 — cos^opt)/2. Fig. [T31 shows that this linear relationship 
for cos 6opt is a valid description for s-channel single top events at the parton level. Moreover, 
in contrast to the helicity basis, the top quark is almost completely polarized in the "optimal 
basis" at parton level, and the 0{as) corrections don't change this picture very much. 

However, in this case the event reconstruction itself has a significant effect on the dis- 
tribution. The reconstructed cos 6opt distribution also shows a drop-off, in this case at high 
cosOopti which is due to the r] cut on the lepton. 

To better quantify the change in polarization, it is useful to define the degree of polar- 
ization T) of the top quark. This is given as the ratio 

V = 

iV„ + iV+ ' 

where (A^+) is , the number of left-hand (right-hand) polarized top quarks in the helicity 
basis. Similarly, in the optimal basis, is the number of top quarks with polarization 
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Figure 15: Top quark polarization in the "optimal basis" using the fuh parton information (left) and 
after event reconstruction (right) with selection cuts, comparing Born level to 0{as) corrections. 

against (along) the direction of the anti-proton three momentum in the top quark rest frame 
ft. The angular distribution is then given by 



da 



1 + cos 6* 



a d{cos6) 



N_ - 

id 



N+ 2 
- D cos 6j) 



+ 



cos 9 



A simple algebra leads to the following identity: 

V = -3 



^ da , 
X — ^ dx , 

1 



adx 



(10) 



where 



da 
adx 



is the normalized differential cross section as a function of the polar angle x. 



Here, x denotes cos 6hei in the helicity basis and cos 6opt in the "optimal basis". Based on 
the degree of polarization V, we can easily get the spin fractions J^± as: 
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N_^+N+ 2 

Note that JF_(JF_,_) is the fraction of left-handed (right-handed) polarized top quarks in the 
helicity basis. Similar, in the optimal basis, JF_(JF_,_) is the fraction of top quarks with 
polarization against (along) the direction of the anti-proton three momentum in the top 
quark rest frame. 

We can also define the asymmetry A of the distribution as 

jli da{cos6) — da{c 



A 



[cost 



COS I 



lo M 



cos( 



(11) 



98 







V 




A 






LO 
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Helicity Basis: 


Parton level 


0.63 


0.48 


0.58 


0.82 


0.74 


0.79 


0.32 


0.24 


0.29 




Reconstructed events 


0.46 


0.17 


0.37 


0.73 


0.58 


0.68 


0.26 


0.11 


0.21 


Optimal basis: 


Parton level 


-0.96 


-0.83 


-0.92 


0.98 


0.92 


0.96 


-0.48 


-0.42 


-0.46 




Reconstructed events 


-0.48 


-0.28 


-0.42 


0.74 


0.64 


0.71 


-0.24 


-0.14 


-0.21 



Table III: Degree of polarization D, polarization fraction and asymmetry A in s-channel single 
top events. Here, T corresponds to T- in the helicity basis for left-handed top quarks and to Tj^ 
in the optimal basis for top quarks with polarization along the direction of the anti-proton three 
momentum. 



It is easy to check that without imposing any kinematical cuts, D = 2A. Furthermore, the 
ratio of top quarks with spin along the basis direction will be = 0.5 — A when no cuts are 
applied. However, when cuts are imposed, the two relations break down. 

Table ITTTl shows that the relationship V = 2A indeed holds at parton level (within round- 
ing errors) and is still approximately true at 0{as). We present our result at the parton 
level, both before and after event reconstruction. We found that at the parton level stage 
the "optimal basis" presents a significant improvement over the the helicity basis at both 
the Born level and NLO. The top quark is almost completely polarized in the optimal basis. 
However, after event reconstruction both bases give almost the same spin fractions. The 
gain from using the optimal basis is lost again and the helicity basis is equivalent to the 
optimal basis in terms of polarization. This is a result of the jet reconstruction and lepton- 
jet separation cut and of the constraint on the W mass when reconstructing the neutrino. 
Furthermore, the 0{as) corrections reduce the polarization already at the parton level in 
both the helicity and optimal basis. 

Before concluding this section, we note that the relative amounts of production- and 
decay-stage emission depend sensitively on the kinematical cuts applied. Also, we do not 
show the transverse momentum distribution of the t b pair or their azimuthal angle separation 



because those distributions are sensitive to multiu 
reliably predicted using a resummed calculation |4 1 . 



e gluon radiations, and they can only be 
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D. Distributions for three-jet event 



As shown in Fig. |3J a large fraction of the events passing the loose selection cuts contains 
three jets. In this section we focus on the properties of the additional jet. 

1. Kinematical Distribution of the Extra Jet 

Initial- and final-state emission of additional gluons occurs before the top quark goes on 
shell and can thus be considered as "production-stage emission" , while decay-stage emission 
occurs only after the top quark goes on shell. In principle, an event with an extra jet can thus 
be classified as production-stage or decay-stage by looking at the invariant mass of the decay 
products. In pro duct ion- stage emission events, the IV-boson and 6-quark momenta will 
combine to give the top momentum. In decay-stage emission events, the gluon momentum 
must also be included to reconstruct the top momentum. This interpretation is exact at 
the parton level in the narrow width approximation. Finite top width effects can blur it 
due to interference between production- and decay-stage emission. This classification is 
nevertheless still useful in our case because the top width of 1.5 GeV is small compared to 
the hard gluon Et cut imposed in the MC calculations. It should be kept in mind that in 
an experiment, the production-decay distinction is further blurred by the experimental jet 
energy resolution and ambiguities associated with properly assigning partons to jets and the 
like. 

Fig. ^1 shows the transverse momentum distribution as well as the pseudo-rapidity dis- 
tribution for the third jet in 3-jet events. This jet corresponds to the gluon in about 70% of 
the events after the loose set of cuts. 

Note that production-stage emission is dominant over decay-stage emission because the 
decay contribution is determined not by the collider energy, but by the phase space of the 
top quark decay. As expected, the Et distribution is steeply falling for all contributions, but 
it extends to much higher Et values for production-stage emission. The smaller values of Et 
to which decay emission is constrained are again a consequence of the top decay kinematics. 
Recall that the top quarks are produced with relatively modest transverse momentum (cf. 
Fig. ITT|) . so that gluons from the decay do not receive much of a boost. Note also that 
an increase in the Et cut on the jet would result in a further reduction in relative size 
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Figure 16: Transverse momentum (left) and pseudo-rapidity of the third jet after selection cuts for 
the various 0{as) contributions. 
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Figure 17: Distance between the third jet and the best-jet (left) and distance between the third 
jet and the non-best-jet (right) after selection cuts for the various 0{as) corrections. 

of the decay contribution compared to production. Fig. also shows the distribution in 
pseudo-rapidity of the extra jet. Initial-state emission is relatively fiat in pseudo-rapidity, 
as compared to the more central decay emission. This is consistent with our intuition that 
decay-stage radiation, being associated with the final-state particles - which tend to appear 
in the central pseudo-rapidity region - is also lilcely to be produced centrally. But this decay 
contribution is small and pro duct ion- stage radiation dominates even in the central region. 
Also, final-state emission and decay emission have a tendency to follow the direction of the 
6-jet and fe-jet respectively, hence the decay 0{as) contribution peaks at a positive t] value 
while the final state 0{as) contribution peaks at a negative t] value, cf. Fig. [7| 

This tendency of decay-stage radiation to be associated with the final-state b quarks 
might lead one to expect that if the extra jet is "near" the 6-jet it should be included in 



the mass reconstruction, and if it is not then it should be excluded. Fig. El which shows 
the distribution in AR between the extra-jet and the best-jet, confirms that the decay-stage 
radiation peaks close to the best-jet, and pro duct ion- stage radiation peaks farther away. 
Unfortunately, the production contribution is so large that it dominates even at the low 
AR cutoff. A higher Et cut on the jet would make this situation even worse. The best 
choice of what is "near" the b quark will therefore balance the competing effects of decay 
emission falling outside the cone and production emission falling inside the cone. In Fig. El 
the equivalent distribution in AR between the extra jet and the non-best-jet is also shown, 
where the pro duct ion- stage radiation peaks close to the non-best-jet due to the effect of 
gluon radiation in the final state. 

Tuning a prescription for dealing with the extra jet in s-channel single top events is further 
complicated because effects of multiple emission, hadronization, and detector resolutions will 
affect the result. 

Finally, we point out that the above conclusion does not strongly depend on the choice of 
jet algorithm. We have checked that at the parton level using the Durham kr algorithm [l^ 
l2l| leads to a similar conclusion on the relative importance of the production- and decay- 
stage gluon emission. 



2. Angular Correlation Between the Extra Jet and the Best-Jet 

While we cannot use the angular distance between the gluon jet and the 6-jet to dis- 
tinguish production-stage emission from decay-stage emission, it is nevertheless possible to 
separate the two using the best-jet algorithm. 

The left-hand side of Fig. El shows the angular correlation cos6 between the gluon and 
the b quark at parton level before cuts. The right-hand side of the same figure shows the 
same correlation after event reconstruction between the third jet and the best-jet. In this 
case there is a clear separation between pro duct ion- stage and decay-stage emission and the 
best-jet algorithm can be used to separate the two because only those events are included 
in the plot for which the best-jet algorithm consists of exactly two jets. 
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Figure 18: Angular correlation cosO between the gluon and the b quark before any selection cuts 
using the full parton information (left) and between the third jet and the best-jet after selection 
cuts (right). The solid line shows all 0{as) contributions except for the decay part, while the 
dashed line shows only the 0{as) decay contribution. 

V. SINGLE TOP PRODUCTION AS BACKGROUND TO SM HIGGS 
SEARCHES 



The s-channel single top quark process also contributes as one of the major backgrounds 
to the SM Higgs searching channel qq WH with H bb. In this case it is particularly 
important to understand how 0{as) corrections will change distributions around the Higgs 
mass region. 

Because of the scalar property of the Higgs boson, its decay products b and b have 
symmetric distributions. Fig. shows the invariant mass distribution of the 6-jet, fe-jet 
pair. For a Higgs signal, this invariant mass of the two reconstructed fe-tagged jets would 
correspond to a plot of the reconstructed Higgs mass. Thus, understanding this invariant 
mass distribution will be important to reach the highest sensitivity for Higgs boson searches 
at the Tevatron. The figure shows that at 0{as), the invariant mass distribution not only 
peaks at lower values than at Born level, it also drops off faster. This change in shape is 
particularly relevant in the region focused on by SM Higgs boson searches of 115 < m^j < 
130 GeV where the Higgs associated production cross section is also at the fb level. In 
particular the SDEC contribution, while small overall, has a sizable effect in this region of 
the invariant mass and will thus have to be considered in order to make reliable background 
predictions for the Higgs boson searches. 
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Figure 19: Invariant mass of the (6-jet, fe-jet) system after selection cuts, comparing Born level to 
0(qs) corrections. 
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Figure 20: Angular distance cos9 between a 6-tagged jet and the (b jet, b jet) system after selection 
cuts, comparing Born level to 0{as) corrections. 



The higher order QCD corrections are evident in the distribution of cos6 for the two 
6-tagged jets, cf. Fig. HUJ Here 9 is the angle between the direction of a 6-tagged jet 
and the direction of the (&-jet, 6-jet) system, in the rest frame of the (6-jet, fc-jet) system. 
Experiments cannot distinguish between the b- and the 6-jets, we therefore include both 
the 6-jet and the 6-jet in the graph. This distribution is generally flat at Born level, with 
a drop-off at high cos6 due to jet clustering effects, and a drop-off at negative cos6 due 
to kinematics. The 0{as) corrections change this distribution significantly and result in a 
more forward peak of the distribution. In particular the SDEC contribution flattens out the 
0{as) corrections and therefore the s-channel distribution resembles more that of the Higgs 
boson signal. This becomes important when the shape of this distribution is used to extract 
the Higgs boson signal from various backgrounds. 



VI. CONCLUSIONS 



We have presented a next-to-leading order study of s-channel single top quark events at 
the Tevatron, including 0{as) QCD corrections to both the production and decay of the 
top quark. To obtain an accurate prediction of the inclusive rate of s-channel single top 
production, a "modified" narrow width approximation was adopted to link the production 
of the top quark with its decay (and thus preserve top quark spin information) instead of 
the usual narrow width approximation. In the former, a Breit-Wigner distribution of top 
quark mass is generated in the phase space generator of the Monte Carlo program, and the 
generated mass is also used to calculate the scattering matrix elements. Hence, the effect of 
the top quark width is included. In the latter, a fixed value of top quark mass is used. The 
impact of kinematical cuts on the acceptances has been studied for several different sets of 
cuts. We found that the difference between Born level and NLO acceptances is about 5% 
for two different sets of Pt and t] cuts. This difference becomes significantly larger when 
changing the jet clustering cone size and jet-lepton separation cut Rcut from 0.5 to 1.0. 
Consequently, a constant K-factor cannot be used to normalize Born-level distributions to 
NLO, especially for large values of Rcut- 

We categorize the 0{as) contributions to the s-channel single top process into three 
gauge invariant sets: the initial state corrections, the final state corrections and the top 
quark decay corrections. Keeping track of the different categories facilitates the comparison 
between event generators and exact NLO predictions. 

The 0{as) corrections are significant in size and contribute over 40% of the inclusive 
cross section at NLO, a large fraction of which in events with three reconstructed jets in the 
final state. They also affect the shape of some of the important kinematical distributions 
that might be used experimentally to separate the s-channel single top signal from the 
various backgrounds. Among the higher order QCD corrections, the initial state correction 
dominates over the corrections from the final state and top decay processes in all inclusive 
single particle kinematical distributions. This implies that in particular the initial state soft 
gluon resummation effects need to be modeled properly in single top event generators. We 
have also found that the 0{as) decay contribution, while small in size, has a significant 
impact on several distributions. 

In order to study top quark properties such as the top quark polarization, we need to 



reconstruct the top quark by combining the reconstructed W boson with one of the jets 
in the event. There is an ambiguity in the choice of jet for this reconstruction due to the 
presence of two b jets in the event (the b quark from the top quark decay and the b quark 
produced with the top quark). We use two different methods to resolve this ambiguity: the 
leading jet (highest pt jet) and the best-jet (giving a Wj invariant mass closest to 178GeV). 
Our study shows that the best-jet algorithm identifies the correct jet in about 80% of the 
events while the leading jet is the correct jet in only 55% of the events. We thus choose the 
best-jet algorithm to explore kinematical correlations in the event. 

The top spin correlation is examined in both the helicity basis and optimal basis. At 
parton level, without any kinematical cuts, both bases show strong correlations, with a 
polarization fraction of 98% in the optimal basis and 82% in the helicity basis. The measured 
polarization as well as the difference in fraction of polarization between the optimal basis and 
the helicity basis is reduced after event reconstruction due to the event selection cuts and 
the top quark reconstruction procedure (imperfect W reconstruction). It is further reduced 
by the 0{as) corrections, resulting in a fraction of polarization of about 70% in both basis 
at NLO. 

As a major irreducible background to W^H searches, the s-channel single top quark 
process needs to be well understood in order to discover the Higgs boson or set a bound on the 
Higgs boson mass. We showed that the 0{as) corrections have an effect on distributions that 
are used in W^H searches, in particular in the Higgs mass region 115 GeV < ttlh < 130 GeV. 
The contribution from the top quark decay process, has a sizable effect in this region of 
the invariant mass distribution of the b and b jets. The decay contribution also affects the 
shape of other kinematical distributions, and reliable background predictions in Higgs boson 
searches will therefor need to include higher order QCD corrections. 

In this work we have studied the s-channel single top process at NLO and shown that there 
are significant differences compared to Born level in the cross section and several kinematical 
distributions that are relevant to select single top events at the Tevatron. A more realistic 
study including the analysis of background processes and the comparison of signal rates and 
distributions predicted by various theory calculations will be presented elsewhere. 
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